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1 Introduction 

By now, jump-diffusion processes and more general Levy processes or Levy 
driven processes are widely used in models of the stock dynamics and term 
structure models, under both historic and risk neutral measure. However, only 
discrete time series of empirical data are available, therefore the problem of 
identification of the jump and diffusion components is a non-trivial task. The 
problem of the identification of the diffusion component was considered in de- 
tail in Ait-Sahalia (2002, 2003); in the present paper, the focus is on the jump 
component. Clearly, the presence of jumps is most prominent in prices of contin- 
gent claims near expiry, especially for out-of-the money options, and therefore 
it is natural to extract information about the jump part of the process by us- 
ing empirical data for out-of-the-money European options. Carr and Wu (2003) 
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showed that in the presence of jumps, the prices of out-of-the-money European 
options near expiry are of order r, where r is the time to expiry, and use this 
observation to test for the presence of jumps. However, the results in Carr and 
Wu (2003) are qualitative rather than quantitative. 

The principal aim of the present paper is to derive simple analytical for- 
mulas which can be used for the identification of the jump part of a Levy 
process (under a risk-neutral measure chosen by the market); the straightfor- 
ward generalizations to mean-reverting processes driven by Levy processes (the 
Ornstein-Uhlenbeck type processes), stochastic volatility models with jumps and 
Levy-driven term structure models will be published elsewhere. The first group 
of results of the present paper are formulas for 6, or time decay, of out-of-the 
money European call and put options on a stock, at expiry. Contrary to the 
Gaussian case, 8 for the former (resp., the latter) is shown to be negative if 
there are positive (resp., negative) jumps. In other words, we prove that near 
expiry, these prices are linear functions of time to expiry (up to relatively small 
errors), and calculate the coefficients. The results are explicit analytical ex- 
pressions in terms of the parameters which characterize the process. By using 
these formulas and empirical data for options near expiry as in Carr and Wu 
(2003), it is possible to identify the jump part of a particular Levy model. After 
that, by using the data for options far from expiry, one can identify the first 
instantaneous moment of the process, and calculate the diffusion coefficient by 
subtracting the second instantaneous moment of the jump part of the process. 
This gives a new identification method for the diffusion component. 

By using the formulas for the 0's for the European put and call, we obtain 
approximations for the American put price and early exercise boundary near 
expiry The pricing of American options in the Gaussian case is well understood 
- see e.g. Musiela and Rutkowski (1997), Karatzas and Shreve (1998), Carr 
and Faguet (1994), Carr (1998), Longstaff and Schwartz (2001), Clement et 
al. (2002), and the bibliography therein. In the infinite horizon case, explicit 
analytical formulas are obtained, and in the finite horizon number of 

efficient numerical methods are developed. Although explicit formulas are not 
available in the finite horizon case, several general results are proved for this 
case as well. Consider the American put with the strike price K and maturity 
T, on a non-dividend paying stock; the riskless rate r > is constant. Let 
H(t) = H(r, K,T;t) be the optimal exercise price of the American put. If the 
stock log-price X t = log St follows the Brownian motion, then the asymptotic 
behavior of H(r, K, T; t) is known both near expiry: 

log(ff(r, K, T; t)/K) ~ —a \J(t -T) ln(T - t), (1.1) 

as t — > T — (see, e.g., Barles et al. (1995) and Lamberton (1995); an asymptotic 
formula for the put price is also available), and for small values of r: for any 
t < T, 

H(r,K,T;t)-->0, as r -> +0 (1.2) 

(see, e.g. Musiela and Rutkowski (1997); (|1.2(l is formulated as follows: if the 
riskless interest rate is zero, then the early exercise of the American put is not 
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optimal). Notice that (|l.lfl implies that 

lim H(r, K, T; t) — K, (1.3) 

and that l|1.3|) and l|1.2[l do not agree well when both the riskless interest rate 
and time to expiry vanish. 

In Levendorskii (2003), a variant of Carr's randomization method was de- 
veloped for wide families of Levy processes, and as a by-product, it was shown 
that for many families of non-Gaussian processes used in empirical studies of 
financial markets, the analogs of 11. 3|) and II. 2|) agree much better. Namely, 

lim H(r,K,T;t) = H T (r,K), (1.4) 

where Ht {r, K) depends on some of parameters which characterize the price 
process, and is smaller than the strike price, K. Moreover, it was proved that 
HT(r,K) vanishes with r, and therefore, 

H(r,K,T;t) -> 0, as r -> 0, uniformlyinte[0,T). (1.5) 

Notice that (|1.5|) looks more natural than (|1.2|) , and may be of some importance 
in today's almost deflationary world. Further, if the difference K — Hr(r,K) 
is sizable, then any numerical procedure or parameters' fitting procedure which 
does not take this fact into account, may produce significant errors near expiry 
and strike. Hence, the non-standard behavior of the early exercise boundary has 
practical implications for the parameters' fitting and development of accurate 
numerical methods . 

However, the proof in Levendorskii (2003) has three disadvantages, which 
stem from the use of the Wiener- Hopf method. First, the non-standard behav- 
ior of the early exercise boundary results from long calculations without clear 
intuition. Next, the approximate formulas for the put price are expressed in 
terms of the factors in the Wiener-Hopf factorization formula, hence, they are 
complicated for practical use, and finally, the method in Levendorskii (2003) is 
not directly applicable to a popular class of Variance Gamma Processes, which 
was used in Carr et al (2002), Madan and Hirsa (2003), and Carr and Hirsa 
(2003) for pricing of American put. The generalization of the method for more 
general classes of processes seems impossible. 

Here we prove (|1.4|) and l|1.5l) by using a much simpler argument. We obtain 
an efficient approximation for the put price near expiry, and formulate the re- 
sult about the non-standard behavior of the early exercise boundary in a more 
meaningful form. Also, we extend the result in Levendorskii (2003) to Vari- 
ance Gamma Processes, and in addition, we find a condition on a process which 
ensures that the behavior of the boundary is even more regular than in the 
Gaussian case: there exists a negative constant h° s.t. 

H(r, K, T; t) — K ~ h° ■ (T — t) < 0, &s t —> T. (1.6) 

We prove that the limit of the boundary at expiry is located where —9 of the 
European call is equal to the riskless rate. (If the —9 of the European call at 
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expiry is less than the riskless rate for all spot prices less than the strike, then 
the limit is the strike.) When the corresponding equation is solved, one can 
use the formulas for 6's of the European call and put and obtain an efficient 
approximation to the put price near expiry. The formulas are fairly simple and 
can be used for parameters' fitting purposes near expiry, both for the European 
options and American put. They also can be used to improve performance of 
numerical methods: when using the backward induction, one can start not at 
expiry, where the value function is most irregular but close to it, and use the 
analytical approximation as the basis for the backward induction. 

We demonstrate how large the deviations from the standard behavior of the 
early exercise boundary and put price can be by using the empirical data for 
the American put documented in Carr et al. (2002). We found that the non- 
standard behavior should be observed for reasonable values of the riskless rate: 
typically, the margin between the strike and boundary is 10% or more, and for 
the spot price in the 5%-interval around the strike, the "jump premium" in the 
put price over the pay-off is 0.5r or more, where r is the time to expiry (with 
the normalization t = 1 for a business year). 

The last aim of the paper is to describe the restrictions on the log-price 
process under an EMM, which should be satisfied if the non-standard behavior 
of the early exercise boundary is observed. As (|1.4I) and 11.6( 1 imply, 

various types of the behavior of H{r, K, T; t) near expiry are possible. However, 
for many families of processes used in empirical studies of financial markets, 
namely, Hyperbolic Processes, Normal Inverse Gaussian processes (NIG), and 
Truncated Levy Processes of the extended Koponen's family (a.k.a. KoBoL 
processes, a.k.a. CGMY model) of order v £ (1,2), the behavior of the early 
exercise boundary is non-standard for all parameters' values. This means that 
if in the empirical studies, the early exercise boundary for put options on a 
particular stock exhibits the standard behavior near expiry, then the use of any 
of the processes above to model the log-price dynamics under a risk-neutral 
measure chosen by the market is inconsistent with the data. Normal Inverse 
Gaussian processes, Hyperbolic Processes and KoBoL processes of order greater 
than 1 are processes with sample paths of infinite variation. Notice that the non- 
standard behavior of the early exercise boundary is proved for much wider classes 
of infinite variation Levy processes as well. For processes with finite variation 
jump component, the prominent examples being Variance Gamma Processes 
and KoBoL processes of order less than 1, both the non-standard behavior (|1.4|l 
and "super-regular" behavior ((1.6(1 of the early exercise boundary are possible, 
depending on parameters of the process; there are also borderline cases, when 
we are able only to prove that l|1.3|) holds. According to a recent empirical study 
in Carr et al. (2002), finite variation processes fit empirical data better. 

Even for a finite variation process, if we fix all the parameters of a process 
except for the drift fj,, let r — ► and change /i accordingly so that the stock 
remains a martingale, then below a certain critical value of r (depending on 
the parameters of the process), the non-standard behavior is observed, and as 
r decreases further, the optimal exercise price goes to zero uniformly over the 
interval [0,T). 
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It is natural to ask how robust the result about the non-standard behavior of 
the early exercise boundary is. In Levendorskii (2003) and the present paper, it 
is assumed that dynamics of log-prices X t = log St follows a non-Gaussian Levy 
process under a risk- neutral measure chosen by the market. This simplest choice 
of a non-Gaussian risk-neutral measure for pricing of the American put was em- 
ployed in Carr et al (2002), Madan and Hirsa (2003), Carr and Hirsa (2003), 
where parameters' fitting and computational issues where addressed. Better fit 
can be achieved with stochastic volatility models with jumps but the method in 
Levendorskii (2003) cannot be generalized for stochastic volatility models. The 
method of the present paper can be applied to stochastic volatility models (and 
term structure models), and under certain parameter restrictions, the margin 
between the early exercise boundary and strike can be derived. Further, in Lev- 
endorskii (2003) and the present paper, the formulas for the margin are derived 
in two different approximate models, and in cases which were considered in the 
both papers, the formulas for the margin are the same. The identical results in 
two approximations allows one to hope that the margin in the continuous time 
model is the same. As a proven theorem, we can only claim that the margin in 
the approximate model is not less than in the continuous time model. 

The rest of the paper is organized as follows. In Section |3 we recall the 
definition of Levy processes, and give examples of several families of Levy pro- 
cesses used in empirical studies of financial markets. In Section [21 0's of out- 
of-the-money European put and call options are calculated, and in Section 
the approximate formulas for the American put price and early exercise bound- 
ary near expiry are derived. In Section we consider several families of Levy 
processes used (and suggested for use) in empirical studies of financial markets, 
and discuss the implications of non-standard behavior of prices and early exer- 
cise boundary for the choice of a family of processes, and parameters' fitting in 
more detail. In Section |BJ we derive properties of the early exercise boundary 
for KoBoL processes, for parameters' values documented in Carr et al. (2002), 
and show that the margin between the strike price and early exercise boundary 
can be significant indeed: more than 10%. In the appendix, the technical proofs 
are given. 



2 Levy processes 

Recall that a Levy process is a process with stationary independent increments 
(for general definitions, see e.g. Sato (1999)). A Levy process may have a 
Gaussian component and/or pure jump component. The latter is characterized 
by the density of jumps, which is called the Levy density. We denote it by F(dx). 
Also, a Levy process can be completely specified by its characteristic exponent, 
ip, definable from the equality E [e 4 ^^] = e~ 1 ^^ (we confine ourselves to 
the one-dimensional case). If Xt is a Levy process with finite variation jump 
component, then the characteristic exponent is given by 

2 r+oo 

m = -^+^ 2 + (l-e^)F(dy), (2.1) 
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where a 1 and /i are the variance and drift coefficient of the Gaussian component, 
and F(dy) satisfies 

/ min{l, \y\}F(dy) < +00. 
Jr\{o} 

Equation (|2.1|) is a special case of the Levy-Khintchine formula; for the general 
case, see e.g. Sato (1999). 

Wide families of jump-diffusion processes used in the theoretical and empir- 
ical studies of financial markets are Levy processes with finite variation jump 
component. 

Example 2.1. Let X be a Levy process with the Levy density 

F(dx) = c+X+e^l^^^dx + c_(-A_)e A - x l (0 , +oo) (x)<£r, (2.2) 
where A + > > A_, and c± > 0. Then 

i>\£) = irt - > —7T + t , 2.3) 

2 A + + i£ A_ + i£ 

where a 2 > and /ieR are the variance and drift of the Gaussian component. 
The -0(0 is holomorphic in the strip 3£ G (A-, A+). In order that E[e Xt ] be fi- 
nite, we need to impose an additional condition A_ < —1. For the sake of brevity, 
we will consider below the exponential jump-diffusions with the characteristic 
exponent of the form l|2.3|l . The generalization to the case of Levy densities 
given by general exponential polynomials on positive and negative half-axis is 
straightforward . 

Example 2.2. Variance Gamma Processes (VGP) have been developed and used 
by Madan and co-authors in a series of papers during 90th (see Madan et al. 
(1998) and the bibliography there). The characteristic exponent of a VGP with- 
out the diffusion component can be represented in the form 

= -ifit + c[ln(A+ + i£) - In A+ + ln(-A_ - if) - In(-A_)], (2.4) 

where A + > > A_, c > and // e R. The condition E[e Xi> ] < +00 imposes an 
additional restriction A_ < — 1. 

Example 2.3. Truncated Levy Processes (TLP) constructed by Koponen (1995) 
were used for modeling in real financial markets in Bouchaud and Potters (1997), 
Cont et al. (1997) and Matacz (2000); a generalization of this family was con- 
structed in Boyarchenko and Levendorskii (1999, 2000) and called the extended 
Koponen family of TLP processes. Later, this generalization was used in Carr 
et al. (2002) under the name CGMY-model. As A.N. Shiryaev remarked, the 
name TLP was misleading, and so starting with Boyarchenko and Levendorskii 
(2002a-2002c), we use the name K0B0L processes. 

The characteristic exponent of a K0B0L process is of the form 

^(0 = + cT{-v)[\ + » - (A + + *0' y + (-A-)' y - (-A_ - iO v ], (2-5) 

where v € (0,2), u ^ l,c > 0,A_ < < A + , and /1 6 R. The condition 
i?[e Xt ] < +00 imposes an additional restriction A- < —1. (In Boyarchenko 
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and Levendorskii (1999, 2000, 2002b), the reader can also find a formula for the 
KoBoL process of order v — 1). 

For v £ (Q, 1), the equation l|2.5|l is obtained from l|2.1|l with the Levy density 
given by 

F{dx) = c\ + e x+x l { „ oofi] (x)\x\- , '- 1 dx + c{-\^)e x - x l { o :+oo) {x)x-' / ~ 1 dx, (2.6) 

and cr = (that is, there is no Gaussian component). In the case v £ (1,2), 
instead of i|2.1[l . the general form of the Levy-Khintchinc formula is needed (see 
Sato (1999) and Boyarchenko and Levendorskii (2002b)). KoBoL processes of 
order v £ (0, 1) are finite-variation processes, and the ones of order v £ (1,2) 
are infinite-variation processes. 

Other examples of infinite variation processes are Hyperbolic Processes (HP), 
Normal Inverse Gaussian processes (NIG) and Normal Tempered Stable (NTS) 
Levy processes of order v £ (1,2); in the general classification of Boyarchenko 
and Levendorskii (1999, 2000, 2002a-2002c), HP and NIG are processes of or- 
der 1. Hyperbolic Processes were constructed and used by Ebcrlein and co- 
authors (see e.g. Eberlein and Keller (1995), Eberlein et al. (1998), Eberlein and 
Raible (1999)); hyperbolic distributions were constructed in Barndorff-Nielsen 
(1977). Normal Inverse Gaussian processes (NIG) were introduced in Barndorff- 
Nielsen (1998) and used to model German stocks in Barndorff-Nielsen and Jiang 
(1998); generalization of the class NIG, namely, the class of Normal Tempered 
Stable (NTS) Levy Processes, was constructed in Barndorff-Nielsen and Leven- 
dorskii (2001) and Barndorff-Nielsen and Shephard (2001). As KoBoL processes, 
NTSLP can be of any order between and 2. 

Example 2.4-. The characteristic exponent of an NTS Levy process is of the form 

HO = -itf + S[(a 2 - (/3 + i0 2 Y 12 - (a 2 - P 2 Y /2 }, (2.7) 

where S > 0, a > \/3\, and v £ (0,2). The condition _B[e Xt ] < +oo imposes an 
additional restriction —a + j3 < — 1. With v = 1, we obtain the characteristic 
exponent of NIG processes. 

Remark 2.1. For the explicit calculations in the next two sections, it is 
important that the characteristic exponent of any of Levy processes consid- 
ered above admit the analytic continuation into the complex plane with cuts 
(— £oo,iA_], [iA+, +zoo), where A_ < —1 and A + > 0, and the analytic con- 
tinuation is defined by the same formula (in Example 2.4, A_ = — a + (3 and 
A + = a — (3, and in Example 2.1, the characteristic exponent is analytic in the 
complex plane with two poles at i\- and «A+). This property is unnecessary 
if we want to obtain only qualitative results. For instance, the existence of a 
non-vanishing margin between the strike and early exercise boundary will be 
proved for more general class of Levy processes with the jump part of infinite 
variation. 
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3 Time decay of out-of-the-money European call 
and put options, at expiry 

To simplify the presentation, we normalize the strike price to 1. Let X t be a 
Levy process of any class considered above. Let C{x, t) and V(x, r) be the price 
of the European call and put option, respectively, at time t = T — r, and the 
spot price St = e x . Denote by C(x),x < 0, and V(x), x > 0, the opposite to the 
9 = 9(x,0) of the out-of-the-money European call and put options, respectively, 
at expiry (the option's is the derivative of the option price with respect to 
time): for x < 0, 

C(x) := lim ^Elll (3.1) 

and for x > 0, 

V M ■- lim ^lll. (3.2) 

If the process is Gaussian, then the prices of out-of-the-money European call, 
C(x, t),x < 0, and put, V(x,t),x > 0, vanish faster than any power of r, as 
r — » +0, hence the limits (|3.1|l and (|3.2|l equal 0. The next theorem and lemma 
show that in the case of non-Gaussian Levy processes which we consider, both 
C(x), x < 0, and V(x),x > 0, exist and do not vanish. The intuition is that when 
the jumps are present, the value of waiting of a positive (respectively, negative) 
movement in the price is larger than in the Gaussian case. 

Theorem 3.1. Let X t be any of Levy processes considered above. Then 

a) for x < 0, the limit IS. 1\) exists, and it is positive; 

b) for x > 0, the limit IS. g|) exists, and it is positive; 

c) the time decay at expiry, of out-of-the money European call, C(x), x < 0, and 
put, V[x),x > 0, depends only on the positive and negative jump parts of the 
process, respectively, but not on the drift and Gaussian component. 

Proof. Notice that the intuition behind the independence of the time decay at 
expiry of the Gaussian component is simple: the movements caused by the latter 
are of order r 1 / 2 , on average, and therefore the main contribution to the price of 
an out-of-the-money option, (almost) at expiry, can come from the jumps in the 
corresponding direction only (for a rigorous proof, see the appendix). We start 
with the simplest case of a pure jump process (Example 2.1 with a = 0). Let 
L be the infinitesimal generator of X. Then the price f(X t ,t) of a contingent 
claim with the terminal pay-off g(Xx) can be represented in the form 

f(x, t) = exp[-r(r - L)]g(x + fir). (3.3) 

For the pure jump process, L acts as follows: 



Lf(x)= / (f(x + y)-f(x))F(dy), 
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and since F(dy) < +00 for processes in Example 1.1, we conclude that L 
is a bounded operator in L^R). Hence, if g is bounded (the case of the put 
option: g{x) = (1 — e x ) + ), we can write (I3.3|) in the form 

00 j 

f(x,t) = J2^(L-ry 9 (x + f xr). (3.4) 

If the put option is out-of-the-money, that is, x > 0, then for small r, g(x+fir) = 
0, hence the zero-order term in i|3.4fl is zero, and 

f(x,t) = t(L - r)g(x + fir) + o(r) 
= 7-ig(a;) +o(r). 

By dividing by r > and passing to the limit, we obtain 

V{x)=Lg{x). 

It remains to calculate Lg{x) at x > 0. We have 

/•+00 

£<7(z) = / ((1 - e x+ «) + - (1 - e x ) + )F{dy) 



(1 - e x+ y)+F(dy) 
(1 - e^c+A+e-^+^dy, 
and by calculating the integral, we obtain 

nx) = C -TJ^, z>0. (3.5) 

The case of the call option requires essentially the same calculations but L should 
be regarded as a bounded operator in a space of continuous functions, which grow 
as e x as x — > +00 (for details, see Boyarchenko and Levendorskii (2002b)), and 
the result is 

C(x) = — — , x<0. (3.6) 

— 1 — A— 

Similarly, C(x) and V(x) for more general exponential jump-diffusions can be 
calculated. Moreover, the proof above gives 

/+00 
(e x +y-l)+F(dy), x<0, (3.7) 
-00 

and 

r+00 

V(x)= (l-e x+v )+F(dy), x > 0, (3.8) 
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provided j^_ oQ F{dy) < +00 and 



+00 



e v F(dy) < +00, (3.9) 



(the latter condition is necessary for the stock to be priced). We conjecture that 
<|3.7|l and (|3.8|l are valid for any Levy process satisfying |3.9p . However, we were 
unable to prove (|3.1[) - l|3.2|l in the full generality, and so in the proof of Theorem 
13. II for the other families of Levy processes listed above, we calculate the limits 
C(x) and V{x) directly, by using the formulas for the characteristic exponents. 
In the case of NTS Levy processes and NIG, set A_ = —a + (3, A + = a + (3. 
Now, for any Xt of KoBoL, NTS, and VGP classes, and for z < A- and z > A+, 
define 

*(z) = i[tp(iz + 0) - tp(iz - 0)]. 

Lemma 3.2. Let X t be a VGP, KoBoL, or NTS Levy process. 
Then for x < 0, 

C{x) = (27T)" 1 / y -^-dz, (3.10) 



z(l + z) 
and for x > 0, 

V(x) = (2n)- 1 / — 1 U z. (3.11) 

Jx+ z(l + z) 

Proof. In the appendix. □ 

In the Gaussian case, ^(z) = 0, V z, and since ^ depends linearly on -0, we 
conclude that ^ depends on the jump component of the process only. Hence, 
C{x) 1 x < 0, and V(x), x > 0, depend on the jump component only. The fact that 
the former depends on the positive jumps only, and the latter on the negative 
ones, follows by inspection of the Levy-Khintchine formula: for z > A + (resp., 
z < A_), 'fy(z) is independent of the density of negative (resp., positive) jumps. 

In the appendix, we will calculate ^(z) for each family of processes, and 
derive from (|3.10(l - (|3.11|l the following formulas: 

1) for KoBoL processes of order v e (0, 2), v ^ 1: 

C(x) = - Cr ^ Sh - e ^ ( r + ^V - < 0; (3-12) 
V{x) = r*- X y-\+r dz , x>0 ; (3.13) 

the integrals converge since A_ < — 1 < < A + and v > 0; 

2) for VGP: 



C{x) = c / — dz, x < 0; (3.14) 

> z{l + z) 



V(x) = c / — -dz, x>0; (3.15) 

z(l + z) y ' 
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the integrals converge since A_ < — 1 < < A + ; 
3) for NTS Levy processes 



C(x) = — sin— / dz, x < 0; (3.16) 

" ' , z{z + 1) v ' 





s 




TIP 






sin 






7T 




T 




6 




TIP 






sin 






7T 




~~2 



-a+B 



- xz [(z-[3) 2 -a 2 \ v l 2 



V(x) = -sin— LV , ' - dz, x>0; (3.17) 

I a+B Z(Z+i) 

with v — 1, formulas for NIG processes obtain. The integrals converge since 
-a + [3 < -1, a + fi > 0, and v > 0. □ 

By direct inspection of formulas (|3.12|) ~ (|3.17l) . we conclude that with x = 0, 
the integrands decay as \z\"^ 2 as z — ► ±oo for processes of order v G (0,2); 
if is a VGP or exponential jump-diffusion, then the integrands decay faster 
than |z| e ~ 2 , for any e > 0. Therefore, with x = 0, the integrals diverge if and 
only if the jump component of X t is a process of order v £ [1,2). Further, 
the integrands are positive monotone continuous functions of x. By using the 
Monotone Convergence Theorem, we deduce 

Theorem 3.3. a) For processes of order v G (0, 1), VGP, and exponential 
jump- diffusion processes in Example 2.1, there exist finite limits 

C(-0) : = lim C(x), (3.18) 

X—* — 

V(+0) : = lim V(x). (3.19) 

b) For NIG, and KoBoL and NTS Levy processes of order v G (1,2), the 
limits C(— 0) and V(+0) exist but are infinite. 

c) C is an increasing continuous function, which maps (— oo, 0) onto 
(0,C(— 0)), andV is a decreasing continuous function, which maps (0,+oo) onto 
fP(+0),0). 

Remark 3.1. a) It is possible to derive formulas for C(x), x < 0, and V(x), x > 
0, for the case of Hyperbolic Processes, and that part b) of Theorem 13.31 holds 
in this case. 

b) Part b) can be proved for more general classes of Levy processes. For 
instance, C(— 0) = +oo can be proved if in a neighborhood of 0, the density of 
positive jumps, F + (dy), admits a lower bound via cy~ 2 dy, where c > 0. Then 
there exists a density F^idy) such that / + °° e y F 2 + (d?/) < +oo and 

F+(dy)>F+(dy)+F+(dy), 

where F^(dy) = cy~ 2 e~ v . Clearly, C(x) cannot increase if we replace F + (dy) 
with F^~ (dy) + F^~ (dy) . The argument in the proof of Theorem 13.11 for the 
jump-diffusion case shows that the contribution of F^idy) to C(x) is bounded 
uniformly in x < 0, and the proof for KoBoL processes shows that F^(dy) gives 
a contribution which is unbounded as x — » —0. 
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To finish this section, we list approximate formulas for European call and put 
options near expiry, which follow from Theorem 13.11 and the put-call pairity. Of 
course, these approximations are too inaccurate in a very small neighborhood of 
the strike if the process is of order v > 1. 

For x > 0: 

V(x,t) ~ tV{x); (3.20) 
C{x,t) ~ e x -1+T{r-V(x)); (3.21) 



for x < 0: 



C(x,t) ~ tC{x); (3.22) 
P(x,t) ~ 1 -e x + T(r-C(x)). (3.23) 



4 Early exercise boundary for the American put 
4.1 Finite variation processes 

Assume that C(— 0) and V(+0) exist. They are measures of the influence of 
positive and negative jumps, respectively, and the EMM condition is naturally 
formulated in terms of a, r and C(— 0) and V(+0). As it was shown above, the 
limits (finite or infinite) exist for the families of processes in Examples 2.1-2.4. 

Lemma 4.1. Let X t be a Levy process with finite variation jump component. 
Then e~ rt e Xt is a local martingale if and only if 

2 

r-H- y=C(-0)-7>(+0). (4-1) 

Proof. Since the Gaussian component o~B t and jump component Y f of a Levy 
process are independent, e~ rt e Xt is a local martingale if and only if for any s 
and t > 0, 




whereupon 

e (r-<r 2 /2- M )t = E ^ ^ 

and 

which is 

/+oo 
Pt (y)( e y - l)dy. 
-oo 

(Here pt denotes the transition kernel of the semigroup {T t }t>o, 

/+oo 
p t (z)f(y + z)dz. 
-OO 
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Divide by t: 

/•oo /*0 

t -i( e (r-«r»/a-/.)t_ 1)= / r^^^ -l)dy- t- 1 p t (y)(l- e y)dy, 



JO J -oo 

and pass to the limit t — > 0; the result is 1)4. □ 
We rewrite (14.11) as 

A* - V(+0) = r - C(-0) - cr 2 /2; (4.2) 

for processes without the Gaussian component, 1|4.2[1 is simpler: 

/j,-V(+0)=r-C(-Q). (4.3) 

The first theorem below shows that if negative jumps are present but their 
influence is not very large, and compensated by a positive drift (which clearly 
decreases the value of waiting), then the early exercise boundary converges to 
the strike price. Moreover, the convergence is more regular than in the Gaussian 
case in the sense that l|1.6f) holds. 

Let ft* (r) be the early exercise boundary. The continuous time model can be 
approximated by a discrete time model with small time steps A = T/n. In the 
discrete time approximation, the option can be exercised at t — T, T — A, T — 
2 A, . . . only. Denote by H :— Ht-&{t, K) the optimal exercise boundary at the 
last moment before the expiry, and set ft — logiJ. One should expect that the 
behavior of h = h(A) for small A is a good proxy for the behavior of the optimal 
exercise boundary near expiry. In particular, 

lim h(A) = lim hJr), (4.4) 

A->+0 t->+0 

if the limit in the LHS exists. The proof of a simplified version 

lim h(A) > lim hJr) (4.5) 

is straightforward: as we make the class of admissible stopping times smaller, 
the price of the American put does not increase, therefore the early exercise 
boundary does not go down. We were unable to prove the equality l|4.4|l but 
the inequality (|4.5I) suffices to make the conclusion that if the early exercise 
boundary in the discrete time approximation exhibits the non-standard behavior 
in the limit A — > +0, then the behavior of the early exercise boundary in the 
continuous time model is also non-standard. 

Theorem 4.2. Let X t be a VGP or KoBoL process of finite variation, without 
Gaussian component, or an exponential jump process. Let 

li - V(+0) > 0. (4.6) 

Then 

lim = -V(+0). (4.7) 

a^o A 
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In other words, if— 9 of the out-of-the-money European put, at the expiry and 
strike, is less than the drift, then the slope of the boundary near expiry equals 
the -0. 

Proof. At St-a = H, a put owner must be indifferent between exercising and 
holding the put, therefore h — h(A) is the solution to the equation 

/+oo 
P A(y)(l-e h+ » A+y )+dy 
-oo 

(the LHS is the value of the put if it is exercised, and the RHS is the expected 
present value of the put if it is kept alive), equivalently, 

,o 

1 - e h = er rA / PA (y-h- nA)(l - e y)dy. (4.8) 

J — oo 

Assume that for small A > 0, h(A) + fiA > 0, then > h(A) > -/iA, hence 
h(A) — » 0. By using the Taylor formula and dividing by A, we obtain 

-A= [ A- x p A (y-h- n&)(l-eV)dy + 0(&). (4.9) 

^ J — oo 

Then we pass to the limit A — > 

lim ~ fe f A ^ = V(+0). 

Condition V(+0) < /i is necessary lest the assumption h(A) + /u,A > lead to a 
contradiction; and the condition (|4.ti|l suffices for the inequality h(A) + [iA > 
to hold for small A. □ 

Now assume that (|4.6t fails. This happens when either the drift is negative or 
the influence of negative jumps is too large. Hence, the value of waiting should 
increase, and the optimal exercise boundary should be lower. As we will see, 
apart from the borderline case \i = V(+0), the value of waiting becomes so large 
that h — h(A) remains bounded away from 0, as A — > 0. 

Theorem 4.3. Let X t be a VGP or KoBoL process of finite variation, without 
Gaussian component, or an exponential jump process. 
Then a) if fi — V(+0), then the limit 

h=mxi Q h{A) (4.10) 

exists, and it is equal to 0; 
b)if 

fi-V(+0)<0, (4.11) 
then the limit exists. It is negative, and can be found from the equation 

r = C(h), (4.12) 

which has a unique solution. 
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In other words, if —9 of the out-of-the-money European put, at the expiry 
and strike, is bigger than the drift, then the limit of the early exercise boundary 
at expiry is less than the strike. It is located where —9 of the out-of-the-money 
call, at expiry, is equal to the riskless rate. 

Proof. We start with part b). First of all, we notice that due to 1)4.3(1 . (|4.11() is 
equivalent to 

C(-0) > r, (4.13) 

therefore on the strength of Theorem l3.3l c), the solution h to (|4.12() exists and it 
is unique. Secondly, the reader may be surprised that the answer is formulated 
in terms of the density of positive jumps, but the proof shows that this is a 
consequence of the put-call pairity (for European options). Since e ~ rt+Xt is a 
(local) martingale under the risk-neutral measure, we have 

/+oo r + oo 

PA(y)e x+ » A+v dy= / PA (y-x~[iA) e ydy = e x+rA . 
-oo J —oo 

Hence, we can rewrite the RHS in ()4.8|) as 

(r+oa r + oo 

/ PA (y -//A - h)(l - e y)dy - / PA (y — fiA — h)(l - e^)dy 
J -oo JO 
r+oo 

= er rA - e h + e- rA / Pa{p — f-A — h){e v - l)dy 
Jo 

(essentially, this is the put-call pairity), and therefore, (|4.8|l can be written as 

f+oo 

e rA -l= / p A (y-fiA-h)(ey-l)dy. (4.14) 



We see that if /iA + h(A) is negative for small A, then h depends only on the 
upper tail of the probability density. By dividing by A and passing to the limit, 
we obtain l|4.12|l . Thus, part b) has been proved. 

In part a), /i = V(+G), hence from l|4.1|l . r = C(— 0), and if we assume that 
for some sequence A n — > 0, the sequence h(A n ) remains bounded away from 0, 
then by dividing by A = A n in 14.14JI and passing to the limit, we obtain an 
inequality, contradiction. □ 

In the following theorem, we allow for a Gaussian component. 

Theorem 4.4. Let X t be a VGP or KoBoL of order v 6 (0, 1) 7 or exponential 
jump- diffusion process, with a Gaussian component. 

Then a) if r > C{— 0), then the limit h exists. It is negative, and can be found 
from $~Wj) ; 

b) if r < C{— 0), then the limit h exists, and it is equal to 0. 

Proof. By repeating the proof of Theorem 14.31 we obtain a) , and similar argu- 
ments show that under condition r < C{— 0), an assumption that the behavior 
of the boundary is non-standard, leads to contradiction. This gives b). □ 
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We see that in the case r < C(— 0), the result is weaker than in the pure 
jump case, where we can single out the case when the boundary is more regular 
than in the Gaussian case. 

4.2 Infinite variation processes 

According to part b) of Theorem 13.31 now C(x) — * +oo as x — » —0. By using 
this observation, and repeating the proof of part a) of Theorem 14.31 we obtain 

Theorem 4.5. Let X t be a NIG or KoBoL or NTSL process of infinite variation. 

Then the limit h = liniA^o h(A) exists. It is negative, and can be found from 
the equation 

Remark J^.l. If the density of positive jumps satisfies the conditions in Re- 
mark 3.1 b), then C(— 0) = +00, and repeating the proof of Theorem 14.31 we 
obtain that h is negative. 

4.3 Optimal exercise price when the riskless rate vanishes 

Let X t be a process from any of families in Examples 2.1-2.4. Let r — > +0 but 
the Levy density and the diffusion coefficient of the Gaussian component remain 
fixed. Then the drift of the process changes with r and converges to a finite 
value. By the direct inspection of 14.12f> . we conclude that if the early exercise 
boundary exhibits the non-Gaussian behavior i|1.4fl near expiry, then the early 
exercise price near expiry vanishes with r as well. The same argument as in the 
Gaussian case shows that for any t > 0, h*{r) = /i*(t, r) — > +0 as r — > +0 
(see Hl-20 ). therefore we conclude that the optimal exercise price tends to zero 
uniformly on the interval [0,T). 

Notice that if X t is a VGP or a K0B0L process of order v G (0,1), or 
exponential jump-diffusion (as in Example 2.1), then for large r, l|4.12|l has no 
solution, but for sufficiently small r, equation 14.12|) has a unique solution, and 
this solution tends to together with r. This means that for sufficiently low 
levels of the riskless rate, all the processes considered above should lead to the 
non-standard behavior of the early exercise boundary. 

4.4 Approximate formulas for the American put price near 
expiry 

By using the approximate formulas for the European put and call near expiry 
at the end of Section 3, we obtain approximate formulas for the American put 
price near expiry: 

( l~e x , x<h 
v(x,t) ~ I l-e x + t(C(x) -r), h < x < (4.15) 
[ tT{x), x>0 

Notice that this approximation is too inaccurate very close to the strike, espe- 
cially for infinite variation processes. 
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5 Implications for parameters' fitting 

The behavior of the early exercise boundary near expiry can be used to de- 
termine the type of a Levy process, which can be used to model the log-price 
dynamics under the risk-neutral measure. This choice of a risk-neutral measure 
was used in Carr et al. (2002), Hirsa and Madan (2003) and Carr and Hirsa 
(2003). Empirical studies of financial markets clearly indicate the presence of 
both downward and upward jumps (under historic measure), therefore if the 
boundary exhibits the standard behavior, then one may not use simultaneously 
a process of order v 6 [1,2) under the historic measure and Levy process under 
a risk-neutral measure (chosen by the market). Indeed, as it is shown in Carr et 
al. (2002), p. 312, the difference between the historic and risk-neutral Levy den- 
sities must be of finite variation, therefore a process of infinite (finite) variation 
under the historic measure must remain a process of infinite (finite) variation 
under an EMM. Processes of order v < 1 can be used for modelling of log-price 
dynamics both under historic and risk-neutral measure but additional restric- 
tions on the density of the upward jumps must be imposed. Another possibility 
is to use one-sided KoBoL processes with densities of positive jumps of an or- 
der different from the one for positive jumps (see Boyarchenko and Levendorskii 
(1999, 2000, 2002b)). The characteristic exponents of these strongly asymmetric 
KoBoL processes are of the form 



where c± > 0, v± < 2, A_ <0< A+, and v± ^ 1 (for the formulas in the case 
v± = 1, see Boyarchenko and Levendorskii (1999, 2000, 2002b)). The analogs of 
<pn2|l and (|3~T3|) are 



If a finite variation process is chosen for modelling of the American put on a 
particular stock, then the choice of parameters of the process is constrained by 
the observed behavior of the early exercise boundary: 

(i) if the non-standard behavior i|1.4|) is observed, then the constraint is r < 



(ii) if the "super-regular" behavior H1.6(l is observed, then the constraint is 
/j > V(+0) (and probably, there is no Gaussian component); however, as 
empirical examples in the next section show, this possibility is unlikely to 
realize; 

(iii) if ft, = but the "super-regular" behavior (|1.6|) is not documented, then 
the constraint is r > C{— 0), and the Gaussian component is admissible. 



i(l!; + C + T(-v + )[\ + »+ - 

+c_r(-i/_)[(-A_) 



(A+ + iHT + ] 

-(-A- -if) 




C(-0); 
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Still another possibility is to assume that under an EMM chosen by the market, 
the rate of decay of the density of positive jumps is very large; then h = for 
processes of order v < 1, and h is negative but very close to for processes 
of order v > 1. For KoBoL, this means that the steepness parameter A_ must 
be large in modulus. If under the historic measure, the A_ is not large, this 
interpretation presumes that the agents in the market are very risk averse. 

The formulas for 9 of out-of-the money European call and put options, at 
expiry, can also be used to infer the type of an appropriate family of processes, 
and parameters of the process. In particular, if 9 at expiry quickly grows as 
the spot price approaches the strike price, then processes with infinite variation 
jump component must be used, otherwise the processes with finite variation 
jump component. 

6 Empirical examples 

The aim of this section is to demonstrate how large the margin between the 
strike and early exercise boundary and the "jump premium" in the American 
option price over the payoff 1 — e x can be for realistic values of parameters. 
We use risk-neutral parameters' values of the KoBoL processes with a diffusion 
component, which were obtained for several stocks in Carr et al. (2002) (Table 
3 on p. 327 in op. cit.). 

Example 5.1. Only in one case: ibmllll, the process is of order v > 1: 
the risk-neutral parameters are v = 1.0102; c = 0.42; A + = 4.37; A_ = —191.20, 
and a = 0.428 (in Carr et al., these parameters are denoted Y,C,G,—M, and 
rj). Since v > 1, the non-standard behavior of the early exercise boundary is 
guaranteed but the gap between the strike price and boundary is small unless 
the interest rate is very small indeed. It is seen from Figure 1, where we plot 
the graph of h as a function of r, that for r = 0.02, h is about 0.006, and for 
r = 0.005, h is close to 0.01, so the gap is clearly seen. In Figure 2, we plot 
the graph of 9 at expiry, for out-of-the-money European call and put (graph of 
C(x), for small x < near zero, and put, V(x), for x > 0). Both are unbounded 
in a neighborhood of 0. Since the density of positive jumps decays very fast 
(the steepness parameter A_ = —191.20 is very large), the time decay of out- 
of-the-money puts is very small unless the spot price is very close to the strike. 
The density of positive jumps decays faster than the one of negative jumps 
(— A_ > A + ), therefore the 9 for the out-of-the-money call is smaller than the 
one for out-of-the-money put, for the same distance from the strike. 

Example 5.2. For amznl014, the parameters are a = 0.0684; v = 0.3072; c = 
4.60; A + = 1.78; A_ = —6.29. The fairly large c indicates that the intensity of 
jumps is large, and relatively small values of A + and A_ imply that there are 
many large jumps, especially negative ones (recall that A + characterizes the rate 
of decay of the density of negative jumps). Therefore, even not very close to 
expiry, the owner of an out-of-the-money European option can expect that with 
a non-negligible probability the option will become in-the-money. Hence, the 
price of these out-of-the-money options is sizable, especially the one of the put. 
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Figure 1: Limit of log of the early exercise boundary at expiry, as a function of 
the riskless rate. Parameters (ibmllll): v = 1.0102; c = 0.42; A + = 4.37; A_ = 
-191.20, a = 0.428. 
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Figure 2: Time decay of the out-of-the-money call and put, at expiry. Parameters 
(ibmllll): v = 1.0102; c = 0.42; A+ = 4.37; A_ = -191.20, a = 0.428. 

Figure 3 demonstrates this effect. We plot graphs of C and V in a neighborhood 
of zero. Further, C(— 0) = 1.7851, therefore for reasonable r, we have r < C(— 0), 
and the behavior of the early exercise boundary is non-standard. Since the 
intensity of large negative jumps is significant, the value of waiting is very large, 
and the gap between the strike and optimal exercise boundary is also large. 
From H4.15fl and Figure 3, we can conclude that even at 5% below the strike, 
the option value of waiting (the difference between the American put price and 
the pay-off 1 — e x ) is of order 0.5r; r = 1 is one year. 

In the next two Figures 4 and 5, we plot the graph of C farther from the 
strike, and the graph of h as the function of r. In Figure 6, we plot both 
the early exercise boundary, at expiry, H = exp(h(r)), as a function of r, and 
that of the early exercise boundary far from expiry, H« = exp(/i*(r)) (the early 
exercise boundary in the infinite time horizon case) . Similar results are valid for 
the other stocks documented in Carr et al. (2002). 

7 Conclusion 

For several families of Levy processes used in empirical studies of financial mar- 
kets, we have proved that the time decay of out-of-the-money European options 
at expiry is not negligible, and derived explicit formulas. These formulas can 
be used to make numerical calculations near expiry more accurate - both for 
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Figure 3: Time decay of out-of-the-money European call and put, at expiry. 
Parameters (amznl014): v = 0.3072; c = 4.60; A+ = 1.78; A_ = -6.29, a = 
0.0684. 
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Figure 4: Time decay of out-of-the-money European call, at expiry. Parameters 
(amznl014): v = 0.3072; c = 4.60; A + = 1.78; A_ = -6.29, a = 0.0684. 
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Figure 5: Limit of log of the early exercise boundary, at expiry. Parameters 
(amznl014): v = 0.3072; c = 4.60; A+ = 1.78; A_ = -6.29, a = 0.0684. 
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Figure 6: Graphs of the early exercise boundary, at expiry (solid line), and 
far from expiry (dotted line), as functions of the riskless rate. Parameters 
(amznl014): v = 0.3072; c = 4.60; A+ = 1.78; A_ = -6.29, a = 0.0684. 
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European options, and American ones. By using the same formulas, we proved 
that the optimal exercise boundary for the American put is separated from the 
strike price by a non- vanishing margin on the interval [0,T), and that as the 
riskless rate vanishes, the optimal exercise price goes to zero uniformly over the 
interval [0,T). The last result may be especially interesting since the current 
levels of interest rates are rather low. Furthermore, we discussed the restric- 
tions on the process under the risk-neutral measure which the non-Gaussian 
behavior of the early exercise price near expiry imposes, and suggested a natural 
asymmetric extension of the KoBoL family (a.k.a. CGMY model), which can 
be used for any type of the behavior of the boundary. We showed that for the 
risk-neutral parameters obtained for several stocks in Carr et al. (2002), the 
behavior should be non-standard, and the margin between the strike and the 
early exercise boundary be quite sizable: more than 10%, up to expiry. The 
jump premium over the payoff is shown to be sizable as well. 



A Technical proofs 
A.l Proof of Lemma 13.21 

First, we consider processes without the Gaussian component, and in the end 
prove that the results remain the same after an addition of a Gaussian compo- 
nent. 

Further, we assume that A_ < — 1; the results for the case A_ = — 1 obtain 
by passing to the limit A_ | — 1 in the formulas for the case A_ < — 1. For any 
<j G (A_, —1), the price of the European call option is given by 

C(x,r) = (2.)- r +t ° ^ r r(r + m)] ^ (AJ) 



— X - in 



-i£(l-*0 



(see e.g. Boyarchenko and Levendorskii (1999, 2002b)). Clearly, the limit C(x) 
is independent of r, therefore we may set r = 0. The characteristic exponent of 
a process without the Gaussian component is of the form 

*/>(£) = + 0(0, 

where </>(£) grows as in the complex plane with the cuts (— ioo,iA_] and 
-Moo), for processes of order v G (0, 2). in the complex plane outside two 
poles £ = i\±. For VGP, grows as In \z\, and for exponential jump processes, 
decays as 
Set x' = x + [it, and write l|A.l(l as 



Let X t be a KoBoL or NTS Levy process of order v G (0, 1). Then we can find 
S > s.t. p := (1 — 5)/v > 1. Since is bounded away from zero on the 
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line Sy£ = a, the integrand in (|A.2|) admits an upper bound via Ct 1 (1 + |£|) 2 . 
It follows that the integral over > t~ p makes a contribution of order 

r+oo 

2CT- 1 / rf 2 di} = 2Ct p ~ 1 = o(l), as r ^ +0. 

Jt~p 

Since |^)| < (7(1 + IfD", we have r0(£) = 0(r 5 ) = o(l) on the set {£ | 3£ = 
Cj < T_p }- Hence, we may write (|A.2fl as 

r-^,r) - (2TT) - f^ j^I^ e^) 

J-oo+ia - «U 

-oo+icr e ix '^ 



— oc +■?, a 



The a; < is fixed, therefore for sufficiently small r > 0, x' is negative and 
bounded away from zero, and therefore, 5ft(zx'£) is negative and bounded away 
from zero, uniformly in t — ► +0 and £ in the half-plane 5£ < a. This means that 
we can use the Cauchy theorem and push the line of the integration in (jA. 3|> 
down: a — > — oo, and obtain that the integral in (|A.3J| is zero. Hence, 

r-'C(x, r) = (2.)- 1 / ) ^ d£ + o(l). (A.5) 

By using essentially the same argument as in Boyarchenko and Levendorskii 
(2002b), where a similar integral (5.33) was transformed to the integral over 
the banks of the cut (— ioo,iA_] in the complex plane, and after the change of 
variables £ = iz — and £ = iz + on the left and right banks, respectively, 
to an integral over (— oo, A_) (equation (5.37) in Boyarchenko and Levendorskii 
(2002b)), we obtain, for x' < 0, 

/X- -x'zjut \ 
z{1 + [ z) W o(l). (A.6) 

As r — * 0, x' — > x, therefore by passing to the limit in (|A.6() . we obtain (|3.10|) . 

Recall that the calculations above were made for a KoBoL or NTS Levy 
process of order v £ (0, 1). If X t is a VGP, then we make essentially the same 
argument by using any p > 0. If X t is a Hyperbolic Process, NIG, or KoBoL or 
NTS Levy process of order v £ (1, 2), then we start with the integration by part 
in (|A.2(1 . by using 

e lx '^di = (ix'y 1 de lx ' i . 

After that we obtain an absolutely converging integral, and can repeat the con- 
structions above with straightforward changes. In the end, we obtain 



■ x C(a;,T) = {2TTX')- 1 



v ; 1 #(z) ' 



z(l + z) \z{\ + z) 



dz + o(l) 
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By integrating by part back, and passing to the limit r — > +0, we obtain i|3.1U[l . 

To finish the proof of i|3.10[l . it remains to show that the result does not 
change when we add a Gaussian component. Let X t be any of the processes 
above, without a Gaussian component, and let X a .t be a Levy process with the 
same "drift" term and jump density as X t , and with the diffusion coefficient a 2 . 
Let p a t be the probability density of X a , t , and PG,a,t be the probability density 
of the Brownian motion with the volatility a 2 . The pc,a,t is the convolution of 
Pt and PG,a,ti therefore for x > 0, 



r+oo 

t^ 1 C{x,t) = I t~V<t,t(2/ - x)(e v - l)dy 
Jo 



-oc / /> + oc 



t PG,<j,r(z)pr(y - x - z)dz ) (e y - l)dy 

> 

+ 00 Z' + OG 



PG,a,r{z) t p T (y - x - z)(e y - l)dydz. 

00 Jo 

Take p £ (1/2, 1). In the region |z| > r p , the Gaussian density is 0(t n ), for any 
JV, therefore we may replace the outer integral with the integral over z < t p , and 
add o(l). But then on the support of the integrand, for a fixed x < 0, y — x — z is 
bounded away from zero, uniformly in r. Define Cj ump (x) as C(x) but assuming 
that the log-price follows the jump part of the process. Then we can represent 
the inner integral in the form 



(x + z) + o(l), 



and therefore, 



t 1 C(x,t)= / PG,<r,r{z)Cj Ump {x + z)dz + o(l). 

J — 00 

Now, define Cj ump (y) to be Cj ump (x/2) for all y > x/2. Then for small r, the 
integral above does not change but now Cj urnp is a continuous bounded function 
on R. By using the super-exponential decay of the Gaussian density, we can 
restore the , and after that pass to the limit and obtain C(x) = Cj ump (x). 
This finishes the proof of (|3.1UI) . 

The proof of is similar, only we start with the formula for the European 

put: 

J-oo+ia - K) 

where a S (0, A+), and push the contour of the integration up. 



A. 2 Proof of dHHD-dUnD, and (JSHD-dSHD for jump pro- 
cesses with non-zero drift 

If Xt is a VGP, K0B0L or NTS Levy process, then it suffices to calculate &(z) 
and substitute the results in (|3.1U|) and Q3.11[) . 



2G 



1. Let X t be a VGP. We have 

*(z) = ra[ln(A+ - z + iff) + ln(- A_ + z - iO) - ln(A+ -z-iff)- ln(-A_ + z + iO)] . 
If z > A + , then the terms with A_ cancel out, and 

ln(A + -z±iO) = \n(z - A+) ± in. 
Hence, &(z) = —c2n. If z < A_, then the terms with A+ cancel out, and 
ln(-A_ + z ± iO) = ln(z + A_) ± i7r. 

Thus, *0) = c2tt. 

2. Let X t be a KoBoL process. In the case z > A + , 

= icr(-z/)[-(A+ -z + ifff - (-A- + z-iO) 1/ 

+ (A+ - z - iO)" + (-A_ + z + iff) v ] 
= icT{-v)[{\+ ~ z - i0) u - {\ + - z + iQ) v ] 
= icY{-v){z-\+yi{e- 1 ™ -e 1 ™) 
= 2cT{-v)sm{-Kv){z-\ + ) v . 

Similarly, in the case z < A_ , 

#(z) = -2cT(-i/) sin(7iV)(-z + A_)". 

3. Let Xt be an NTS Levy process. If z > A+ = a + j3, then 

*(z) = z<5[(a 2 -(/3 + i(zz + 0) 2 )" /2 -(a 2 -(/3 + i( 2 z-0) 2 ) ,y/2 ] 

= z<S[(a 2 -{(3-z + iff) 2 )"' 2 -{a 2 -(/3-z- z0) 2 )" /2 ] 

= iS(a -f3 + z) u / 2 [{a + f3 - z + iff) v ' 2 - {a + [3 - z - iff) v / 2 ] 

= iS(a -p + zy/ 2 {z -a- fiyl 2 {e 1 ™! 2 - e -™l 2 ) 

= -2^({z-P) 2 -a 2 y/ 2 . 

If z < A_ = —a + 0, then similarly, we obtain the same expression but without 
the minus sign. 

Finally, if Xt is an exponential jump-diffusion process, then the same proof 
as in Lemma 13.21 shows that for the call, it suffices to consider the pure jump 
process, and that in this case, (|A. 5|) holds. Since 



the integrand in i|A.5|) has the only pole £ = iA_ in the plane Q£ < a. We 
push down the line of integration, and on crossing the pole, apply the residue 
theorem. The result is 

r- x C{x,r) = -(2m)- 1 — ^ dg + Q (l) 



+ o(l), 
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and by passing to the limit, we obtain 1|3.6|) . (|3.5|l is proved similarly. 
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